ON THE CONTINUITY OF MULTIVARIATE LAGRANGE 
INTERPOLATION AT CHUNG- YAO LATTICES 



JEAN-PAUL CALVI AND PHUNG VAN MANH 

Abstract. We give a natural geometric condition that ensures that sequences 
of Chung- Yao interpolation polynomials (of fixed degree) of sufficiently dif- 
ferentiable functions converge to a Taylor polynomial. 



1. Introduction 

1.1. Stating the problem. When d+l points ao, . . . , in M converge to a limit 
point a, the corresponding Lagrange interpolation polynomial L[flo, • • • , of 
a function / at the a, 's tends to the Taylor polynomial of / at a to the order d and 
this under the sole assumption that f is d times continuously differentiable on a 
neighborhood of the limit point. This classical result is an easy consequence of 
Newton's formula for Lagrange interpolation and of the mean value theorem for 
divided differences. In this paper, we study a multivariate analogue of this prob- 
lem. We suppose that the points of a multivariate interpolation lattice A of degree 
d in MJ^ converge to a limit point a E and ask under what conditions we can 
assert that the corresponding multivariate Lagrange interpolation polynomials of 
a function / converge to the Taylor polynomial of / at a to the order d ? The 
question is answered for a particular but important class of interpolation lattices, 
the so-called Chung- Yao lattices, see below. 

1.2. A known criterion. In the multivariate case, a simple clear-cut answer can- 
not be expected. This perhaps may be regarded as another consequence of the 
absence of a multivariate mean value equality. We recall a rather general criterion 
(which actually works for hermitian interpolations) which can be found in [QJ. 
Let us mention that the first results which appeared in the literature concerned 
the case (of practical importance in finite elements theory) for which the lattices 
are of the form A*^'^) = U^'\A) where U^'^ is a sequence of linear transformations 
whose norms tend to and A is a fixed lattice. We refer to [|T|| for details and 
references to earlier works. 

We denote by ^'^{R^) the space of polynomials in real variables of degree 
at most d, X" is the monomial function corresponding to the A^-index a, that is 
X"{x) =x"'-- ■x"'' for X = (xi, . . . ,xn) e M^. The length of a is the degree of 
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X", |a| = Y!iLi OCi. We denote by the dimension of the vector space ^'^(M^). 
We have = {^^^) ■ In the whole paper, N >2. 

Theorem 1.1 (Bloom and Calvi). LetA^^^ be a sequence of interpolation lattices 
of degree d in M^. If the following condition holds 

|a| = j+1 ^ limL[A(');X"] =0, (1.1) 

then, for every function f of class C"^'^^^ in a neighborhood of the origin 0, we 
have 

limL[AW;/]=T^(/), (1.2) 

where L[a(^) ; ■] (resp. Tq(-)) denotes the Lagrange interpolation projector at 
the points ofA^'^^ (resp. the Taylor projector at of order d). 

Unfortunately condition ( |1.1[ ) is not easy to verify, especially if the degree of 
interpolation is not small, and it seems difficult to check it on general classes 
of interpolation lattices. Besides, theorem [TT| requires a high order of smooth- 
ness. We point out, however, that although whether the level of differentiability 
required in theorem [TTT] is optimal is not known (in the case of Lagrange inter- 
polation), examples do exist for which convergence does not hold for function of 
class C^^^ but holds for function of higher smoothness, see [jll example 5.4]. 

The aim of this paper is to give a natural geometric condition in the case where 
the interpolation lattices are Chung-Yao lattices. From an algebraic point of 
view, they can be regarded as the simplest interpolation lattices : every point 
is situated at the intersection of hyperplanes chosen among a minimal family 
and the corresponding Lagrange fundamental polynomials are products of affine 
forms. The definition and main properties of Chung-Yao lattices are collected in 
section |2] Our criterion is given and commented in section [3j The proof is quite 
technical and is postponed to the next section. It relies on a remainder formula 
due to Carl de Boor. 

We need very few facts from general interpolation theory. They are recalled 
in the following subsection. 

1.3. Basic facts on interpolation. Let £ be a m-dimensional space of functions 
on and A = {ai ,...,«,„} C M^. We say that A is an interpolation lattice for E 
if for every function / defined on A there exists a unique LeE such that L = / on 
A. Given a basis f = (/i, . . . ,/m) of E, we define the Vandermonde determinant 
VDM(f; A) by 

VDM(f;A):=det(/;-(a,))".^i. (1.3) 
Then A is an interpolation lattice if and only if 



VDM(f;A) 7^0. 



(1.4) 



CONTINUITY PROPERTY OF LAGRANGE INTERPOLATION 3 

Of course, the condition is independent from the choice of the basis f. When 



( |1.4[ ) is satisfied, we have 

m 

L=£/(a,-)l(A,a,v), (1.5) 

where l(A,a,-, •) is the unique element of E which vanishes on A \ {at} and takes 
the value 1 at a,-, 

^ ^ VDM(f;A) - - 

(1.6) 

In the case where E = ^'^(M^) we write L = L[A;/] and call it the Lagrange 
interpolation of / at A. We say that A is an interpolation lattice of degree d. 
The only other case that we consider in this paper is £ = the space 

of homogeneous polynomials of degree d in N variables whose dimension is 

2. Chung- Yao lattices 

We recall the construction of the lattices and of some objects attached to them. 
Despite their apparent simplicity, it seems that these configurations were first 
considered in 1977's Chung and Yao's paper [2]. Here, we essentially follow the 
presentation and notational conventions of de Boor [Q. 

We work in endowed with its canonical euclidean structure. The corre- 
sponding scalar product is denoted by (., ■). 

A set of A'^ hyperplanes H = {ii, . . . , in is said to be in general position 
if the intersection of the A'' hyperplanes is a singleton, that is 

N 



i=l 

If — {x E : {nt^x) = c/}, / = 1, . . . , A^, then H is in general position if and 
only if det(ni , . . . , n^) 7^ 0. 

Definition 2.1. A collection EI of (at least A^) distinct hyperplanes in is said 
to be in general position if 

(1) Every H E (^) — i.e. every subset oiN hyperplanes in EI — is in general 
position (as defined above). 

(2) The map 

HE(^)^^H:=f]ieR'' (2.1) 

is one-to-one. Here and in the sequel we confuse the singleton CleeH^ 
with its element. 



This definition stands at the basis of the following result. 
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Theorem 2.2 (Chung and Yao [2]). Let W be a set of d > N hyperplanes in 
general position in M^. The lattice 



is an interpolation lattice of degree d — N. Moreover, if I & W is given by i = 
{x G : {nf^x) = q} then we have the interpolation formula 

l[0h;/]W= I fm nP4^- (2-3) 

The lattice ©e is called a Chung- Yao lattice (of degree d — N) and the interpo- 
lation formula is called the Chung- Yao interpolation formula corresponding to 
H. In particular, we have 

M^,i>H,x) = Y\P4^^ He(=). (2.4) 

l^%{ne,T^H)-ce \NJ 

When the set of hyperplanes we use is clear, we write instead of ©jj- Of 
course, in ( 2.3| ), different equations for the hyperplanes yield a same formula. In 



the particular case N = \, every set of interpolation nodes may be regarded as a 
(trivial) Chung- Yao lattice. 



As shown by ( |2.3[ ), interpolation polynomials at Chung- Yao lattices are easy 
to compute. Some difficulties, however, must be pointed out. In constructing a 
Chung- Yao lattice, we start from a family of hyperplanes and compute the inter- 
polation points by solving, in principle, m^/ linear systems (of order A''). Besides, 
it is a difficult problem, even in the case N = 2, to decide how to choose the 
hyperplanes if a special requirement is made on the location of the interpolation 
points. For instance, we currently do not know what kind of limiting distribu- 
tion we can obtain with a growing number of Chung- Yao points. We mention 
that an interesting Chung- Yao lattice was constructed by Sauer and Xu [|5l on 
bidimensional disks. 



3. Chung- Yao lattices oe points converging to the origin 

3.1. The convergence theorem. From now on, we shall confuse an hyperplane 
£ with the affine form i{x) = (n,x) — c which defines it, where n is normalized 
so that ||n|| = 1. This abuse of language (each hyperplane has two normalized 
equations) should not create confusion. Boldfaced n will be kept for normalized 
vectors and vectors derived from them. 

Supposing that the points of a sequence of Chung- Yao lattices of same 
degree converge to the origin (or to any other fixed point), we study under what 
conditions the corresponding interpolation operator converges to the Taylor pro- 
jector at the origin. Our main result is summarized in the following theorem. 
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Theorem 3.1. Let d > N. Let s be a sequence of Chung -Yao lattices 
of degree d~N in R^. We assume that Q^"^ is the lattice given by the family of 
hyperplanes 



-1 ' 



with £ 



(n«,.> 



n 



1, / = 1, . . . , J. 

(3.1) 

Consider the following two conditions. 
(CI) All the points of the lattice tend to Q as s ^ °°, that is max{||t?|| : •& G 

(C2) The volumes 



vol(nW,...,nf^ 



1 < /i < Z2 < ■ ■ ■ < < d., 



(3.2) 
are bounded 



of the parallelotope spanned by the vectors n- 
from below, away from 0, uniformly in s. 

If conditions (CI) and (C2) are satisfied then, for every function f of class 
fjd-N+i ^ neighborhood of the origin, we have 



limL[0(-^^/]=T^-^^(/) 



:d-N/ 



(3.3) 



Of course, (3.3 1 holds in every normed vector space topology of ^ ( 

3.2. On condition (C2). The condition on the volume of the parallelotope s is 
equivalent to the following, 



liminf min 



det nr^...,n^^^ 



>0. 



In we have 



where a^j^ ^]0,k[ is the hne angle between the lines 



vo^n-'^^n^''^) = sin(a 



and 



J- 



(3.4) 
(3.5) 

Thus EI satis- 



fies condition (C2) if and only of the angles between any two (distinct) lines in 
remain uniformly bounded from below by a positive constant. An example 
of Chung- Yao lattice of degree 2 in M? and the various parameters involved in 
theorem [3TT are shown in figure [T] 

Conditions (CI) and (C2) are obviously independent. However, when we 
know that the second one holds true, the first one is easily checked as shown 
by the following lemma. 

Lemma 3.2. If(C2) is satisfied then (CI) is equivalent to 



(C3) lim5^oomax,=i , ,^ |c-'*^| = where cf is defined in (li.il). 

Proof We show that (CI) implies (C2). Consider //^ G (^^'') with 6'^ e H^'\ 
From (n-'^-*, — c-'^-* = 0, we get 



< lln 



It? 



Hi' 



//(■') 1 



-s-O, 5-)-oo. 
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dna.} vol]ni,n3) 




Figure 1 . A bidimensional Chung- Yao lattice. 



To show the converse, we observe that if H^''^ - 
nates (xk) of are solutions of the linear system 



(s) 

' I then the coordi- 



N 



Lis) 



J=l, 



and the claim follows from Cramer's formula in which, thanks to condition (C2), 
the denominator remains away from whereas the numerator tends to 0. □ 

3.3. Affine transformations of Chung- Yao lattices. Let ^{x) = L{x) +bhQ 
an affine transformation (isomorphism) of (with L its linear part). If EI is in 
general position so is =Sf (H) := {^(£,) : i= 1, . . . , J} and ^ induces a one-to- 
one correspondence between (^) and {^^^) ■ Moreover if // G {^) then 



^j^^h)='^{^h) and 



In the following theorem we translate the conditions of theorem 3.1 when the 



points of a Chung- Yao lattice are sent to the origin by applying a sequence of 
affine transformations. 

Theorem 3.3. Let EI = {£i, . . . be a fixed collection of d hyperplanes in 
general position in M^, d > N, with, as above, ii = {x G M.^ : (n,,x) — c/ = 0}, 
||n;|| = 1. Let = Ls + bs, 5 G N, be a sequence of affine transformations of 
R^. We set 

e(^) = ^,(e), 5gn. (3.6) 

We consider the sequence of Chung -Yao lattices 0^^-' induced by M.^^\ The fol- 
lowing assertions are equivalent. 

(1) 0'^'') satisfies conditions (CI) and ( C2). 
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(2) There exists a positive constant A such that 

N 

|detL,|-]^||L7^(ni.)|| < A, I <ii < ■ ■ ■ < i^ <d, sen, (3.7) 
;=i 

and 

1 

--h-,d ||L7^(ii,:) 
where Lj^ denotes the transpose of the inverse ofLs. 

Proof. It follows from the normalized equation of J^si^i) together with lemma 



max ^^—jr——-\ci+{ni,L/{bs))\^0, 5 ^ °o, (3.8) 



3.2 Indeed, with £/(x) = {ni,x) — Ci, we have 

^{Q = {x G : {n,,^-\x))-Ci = 0}. 
Since for jc G =Sf(£/), 

0={jVi,L-\x-bs)))-c, = iji,,L-\x))-{c, + {n,,L-\bs))) 

= (L7^(ii,),x)- (c, + (n„L7i(Z7,))), (3.9) 
a normalized equation of =^ (£,) is given by 

^''^"'■^ ^ ^ -{c, + (n„L7i(Z>.))}. □ 



|L7Mn,)lr / llA7Mn/ 



3.4. Examples. In any interpolation lattice of degree 1 is a Chung- Yao lat- 
tice (based on the three distinct lines defined by the interpolation points). More- 
over, any such lattice is the image under an affine isomorphism of the lattice 
:= {(0,0), (1,0), (0, 1)} constructed with the lines of equations (.i{xi,X2) =x\, 
£2(^1,^2) =^2 and £3(^1, X2) =xi +X2 - 1. 

Consider the affine transformations defined by 

where a is a function of t such that lim/^ow(0 = 1^ and the lattice 

0(^) =^s{®) = {{t.t),{t^ + t,t),{t, -t\ + } , t = \/s. 

It is not difficult to see that ©('^ satisfies condition (CI) and (C2). For (C2) we 
use (3.5) and observe that one of the angle is equal to n/l while, thanks to the 
assu mptio n on u, the other two tend to ;r/4 as ? — )• 0. Hence, according to the- 
orem 



3.1 



the corresponding Lagrange interpolation polynomials at &^^> of any 
twice continuously differentiable function / on a neighborhood of converge to 
the Taylor polynomial of /. This example shows that the assumptions of theorem 



3.1[ even in the simple case of theorem [331 are weaker that those given in [V, 



proposition 2.1]. Indeed the assumption 



^ 0, t^O, 



8 



JEAN-PAUL CALVI AND PHUNG VAN MANH 



is required in that proposition whereas it clearly does not hold here since, as is 
easily checked, 

/ I \ , . \ xi — ux\-\- it'^ + t) u — t 

We now give an example showing that convergence to the Taylor projector no 
longer holds, in general, when condition (C2) is not satisfied. We use a compu- 
tation done in [[B example 1.2.]. We fix £ > and define 

0eM = {(O,O),(f,?2+£),(2r,O)}cM2, t=\/s, seN*. 

This lattice satisfies (CI) but not (C2) and it is readily checked that 



2tx\ . 



which clearly does not converge to Tq(X'^2>o)'j = 0as5=l/?— t-oo. The case e = 
shows that that the Lagrange polynomials may converge to a limit different from 
the Taylor polynomial. 

4. Further properties of Chung- Yao lattices and proof of 

theorem 13. ii 

4.1. de Boor's identity. In the following EI always denotes a set of J > hy- 
perplanes in general position in and = 0h the corresponding Chung- Yao 
lattice. We will always assume that 

m = {iu...M- (4.1) 

The elements of EI are ordered according to the indexes. Every subset of EI is 
endowed with the induced ordering. 

If ^ is a subset of — 1 elements in EI, that is ^ G {j^i) , then Hi^x^ is a line 
in which contains d — N +1 points of 0. Indeed, it passes through every 'd-fj 
such that He {^),K CH. The set of these d-N+l points is denoted by 0a:, 

0^ = 0n(^n^j, ^e(^^i)- (4-2) 

Assume that K = {ii^ , . . . , j } with ii < 12 < ■ ■ ■ < in-i- Since the map 

veM^h^det(v,n,,,...,n,-^_;) (4.3) 

is a linear form, there exists a vector, which we denote by n^r, such that 

det(v, n,-; , . . . , J = (v, n/f ) , v G M^. (4.4) 

As defined, the value of n^: depends on the ordering of the hyperplanes of K. 
A different ordering may change in — n^:. It is to avoid further discussion of 
this detail that we assumed we start with a particular ordering of EI and agreed 
that every subset of EI is endowed with the induced ordering. 

Lemma 4.1. If K = . . . ,^!;v_i} then the direction of the line Hi^k^ is given 
by the (nonzero) vector iik- We have WhkW < 1- 
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Proof. The first claim is a consequence of the equations 

=0, j = 1,...,A^-1, 



which follows readily from (4.4). Next, by Hadamard's inequality, the norm of 



the linear form (4.3) is smaller than the product of the ||n,J|'s which is smaller 



than one. Hence, in view of ( |4.4| ), so is the norm of xik- □ 

The vectors n/f play a fundamental role in our proof of theorem [3TT| 
Note in particular that if H E (^) and i E H then we may speak of n^Y- From 
now on, we use H \ ii for H \ {£,}. 

Lemma 4.2 (de Boor's identity). IfH E (JJ) then we have 

^=^«+Ll7^nM^ (4.5) 

where £ denotes the linear part ofi (thus i{x) = (n,jc) if i{x) = (n,x) —c). In 
particular, for every H, the vectors n^y, i E H, form a basis ofM.^. 

Proof See Ep. 37]. □ 

4.2. de Boor's remainder formula. We now recall the definition of multivari- 
ate divided differences. Let Q. be an open convex set in M^, to every set A = 
{ao, . . . , as} C Q. (the points are not necessarily distinct) and / E C'^(f^), we as- 
sociate a 5-linear form on (M/^y defined by 

(M^)^9(vi,...,v.,)^ 

[ao,...,a>i,...,v,]/:= Jd,^...D,J = y'/W(-)(vi, . . . , v,), (4.6) 

[A] [A] 

where /^^^ denotes the s-th total derivative of /, 

and A^. is the standard simplex {t, = (i^i, . . . , i^^) : > 0, YJi=i^i < !}• This 
symmetric 5-linear form is called the multivariate divided difference of / at A. 
Note that, when / E C* (f2) is fixed, the function 

Qf^^ X (M^)' 3 (ao,...,a„vi,...,v,) i — > [ao, . . . | vi, . . . , v,]/ 

is continuous (as a function of its two groups of variables). 
We now state a beautiful error formula due to Carl de Boor. 

Theorem 4.3 (de Boor's remainder formula). Let EI = . . . be a collec- 
tion of d > N hyperplanes in general position in and & = 0e the corre- 
sponding Chung-Yao lattice. For K E (J^j). define the polynomial Pk of 
degree d — N-\-\ by the relation 



10 



JEAN-PAUL CALVI AND PHUNG VAN MANH 



where, as above, I is used for the linear part of I. 

The error between a function f of class C'^^'^+i on a convex neighborhood 
Q. of @ and the Lagrange interpolation polynomial of f at @ is given by the 
following formula. 



f{x)=U@-f]{x)+ £ Pk{x) 



/, xeQ.. (4.8) 



d-N+l 



Recall that for K E (J^J, ®k is the subset formed by the d — N +1 points of® 
lying on the line Hfj^K^ see (4.2). 

Proof See E theorem 3.1.]. □ 

4.3. Some a lgeb raic identities. We now prove two auxiliary lemmas. The first 
one (lemma 4.4) shows that the n/f's, K E (J^J, themselves form a certain 
interpolation lattice. The second one (lemma \4.6^ is a somewhat mysterious 
representation formula for symmetric multilinear forms. 

Lemma 4.4. Let M = {£i, . . . ,£d} be a collection of d hyperplanes in general 
position in with d >N. The set 

is an interpolation lattice for the space J^f^^^^^{MJ^) of homogeneous polyno- 
mials of degree d —N +1. 

Proof. It suffices to prove the following two assertions. 

(1) The cardinality of Y is equal to the dimension of ^^'^^+^(]R^) which 

(rf-v+i) ^ (v-i)- 

(2) For every yik in f there exists Hk E (M^) such that HK{nK) = 1 
but Hk vanishes on T \ {n^:} 

To verify the first point, we just need to check that if K,K' E (J^ J and K ^ 
K' then Hk 7^ n^i. But, if K ^ K' there exists £ eK\K' with i{x) = (ii,jc) - 
c. Assume that K' = . . . j}. Since iUK' is a set of hyperplanes 
in general position, we have det(n,n,j , . . . ,n/^ J 7^ hence, in view of (4.4), 
{n.nic') 7^ 0. On the other hand, since i E K, {n^nx) = 0. Hence iik 7^ njc'. 

As for the second point, for K E (fj^i) , we set 

HK:=PKix)= n #\- (4.10) 

This clearly defines a homogeneous polynomial of degree d — N +1 in sat- 
isfying HK{n.K) = 1. Moreover, if K' E (^^J, K' 7^ K, then we can find I in 
(EI\ /T) n^'. Since t^K, the factor 1{j\.k<^ appears in Hxinic'). However, since 
i E K', i{nK>) = (n,nr) = 0. Hence HxinK') =0. □ 

The interpolation formula corre spond ing to the interpolation lattice — and 
using the polynomials Hjc = in (4.10) — yields the following identity. 
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Corollary 4.5. With the assumptions of the lemma, for every symmetric {d — N + 
I) -linear form ^ on M^, we have 

Hv'-''^')= I PKiv)-H<-''^'). VGM^ (4.11) 
where we use u'^^^^^ := (m, . . . , a) (d — N-\-l times). 

Lemma 4.6. Let M. = {li,. . . ,1^} be a collection ofhyperplanes in in general 
position with d>N. We set 

e,- = {£i,...,£,}, l<i<d, (4.12) 

and 

Then for every symmetric {d —N + \ )-linear form on M^, we /zave 

I /|-^'(^)-<^'(^'"',^/^u£,-,nir''), ^GM''. (4.14) 

In the above formula, we agree that when d — i (resp. / — A^) is not positive then 
X (resp. n^) does not appear in i^-^u^,? w^^)' ^^d when i = d+ \ then x 

and '&k\ju do not appear. Likewise, if the product in the definition of is 
empty then its value is taken to be 1. 

We need the following simple observation. 

Lemma 4.7. Let EI = . . . + 1} be a collection of hyperplanes in in 
general position with d > N. As above, we write = {£i, . . . ,£d}- Let K' G 

{n-i)- VK e andK' ^ K then p|^l(nruwJ = « ^^^'"^ 

5(^1 



emd\K ^("^) 

Proof Take £/ E K' n {Mj \ K) . The fact that G K' gives 

and since ^ i^, £/(n^/u^^+i ) is a factor of (nrui'^+i ) ■ □ 



Proof of lemma 4.6 We prove identity (4.14) by induction on d>N. 
(A) We start with the case J = A^. In that case ( |4.14[ ) reduces to 

H^)= I p^t'\^)H^Kuf,) + I /f'W<^>K) (4.15) 



N-l 

N 



= H^.) + L j-^^Mnu^\iJ. (4.16) 
Since ^ is a linear form, the claim follows from de Boor's identity (|4.5[). 
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(B) We assume that (|4.14[) holds true for d and prove it for J + 1 . Take ^ a 



symmetric (J + 2 — A'^)-linear form. Fix y G and define (^y on (]R^)'^+^ ^ by 
0y(vi,...,vj+i_A,) = ^/)(vi,...,Vrf+i_A,,j). Thus is a symmetric (J+ 1 -A^)- 
Unear form to which we may apply the induction hypothesis to get 



^+1 ,1 . 

I I pl-'\x)^yU-\^Kuuy^'' 



(4.17) 



Putting y = x'm the above expression, we obtain 



(4.18) 



On the other hand, we need to prove 



d-N+2\ 



Expressions (4.18) and (4.19) differ only for / = J-f 1 and i = d + 2. Thus, to 



establish (4.19), it suffices to prove that the term corresponding to J + 1 in (4.18 ) 



equals the sum of the terms corresponding to J + 1 and J -f 2 in (4.19 ), that is 



, TUT ^ ^ / 



d+2-N 



(4.20) 



ForZe (^^J andxGM^, using de Boor's identity (4.5 ) with H = iSTU^^+i, we 
may write 



4+1 w 



-n 



n 



+1 
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Substituting x with the above expression in the left hand side of (4.20), we 
arrive to 

+ I I^^'W /m '^^^ . -^{-r-'^n^mu,,^^). (4.21) 



Now, for K e (^j) , we have 



P^^\x)S^=pI^+\x). 



Hence, the second term on the right hand side of (4.21 1 is 



4+2-N 



(4.22) 



Thus, since K E (5Jl+'), id+i e K means K = K'[J{id+i} with K' E (^^2)' t° 
prove ( |4.20p , it remains to establish 



+1 



We first concentrate on the term Pji^^ (x) j-. — — r on the left hand side of 



(4.23 1. Since i E K v^e have 



p^^hx) 



( 



n 



h{x) 



n 



(4.24) 



(4.25) 
(4.26) 
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Using this expression in the left hand side of (4.23 1, we arrive at 



n 



- E L ^|^kn^,u..,J-<^>(nr^+\n^.u..,,). (4-27) 



{K\e)u£d 



Now, for a fixed K' e (/_^2)' 

using lemma 4.7 for the first equality (we add 
null terms) and corollary C5]for the second one, we get 



Using (4.28 ) in the last term of ( |4.27 1, we finally arrive at 



iix) 



d-N+l 



n 



+1' 



d-N+2 



(4.29) 



which is (4.23 1. This completes the proof of the lemma. □ 

Corollary 4.8. Let M. = {£i, . . . ,id} be a collection of d > N hyperplanes in 
general position in M^. For every function f of class c^^^+i on a convex neigh- 
borhood Q. of the origin in we have 

f{x)-Ti-^{f){x) 

0,...,0,x 



-"5 5 

fl'-A'+l 



Proof. The remainder formula for Taylor polynomial (as a special case of Kergin 
interpolation, see e.g. flU theorem 3]) gives us, 

f{x)-%-\f){x) = [Q^,x\x^]f= J /('^-^+i)(.)(x,...,x). 

d-A'+l d-N+l [0,...,0,x] 

(4.30) 



The corollary then follows directly from lemma 4.6 since, for every a E Q, 
y(t/-A'+i)^^-j ^ symmetric {d — N+ 1 ) -linear form on M^. □ 
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4.4. Proof of theorem 3.1 Let Q.he a neighborhood of the origin on which / 
is of class C^^^ We may assume that 

(i) contains 5(0,i?), the closed euclidean ball of center the origin and 
radius R and, in view of condition (CI), 

(ii) all the points of 0^") = 0jj(.o he in 5(0, i?), 5 G N. 
We set 

M= max ||/('^-^+i)(a)|| <oo, (4.31) 

aeB{0,R) 

where || ■ || here denotes the usual norm of a multilinear form. We use condition 
(C2) in the form given by (|3.4[) taking (|4.4[) into account as follows. 



(iii) There exists 5 > such that 
|(n,-,n^)|>5, Ke 



A^-l 



(4.32) 



(A) We first derive an estimate on the polynomials Pk = defined in (4.1). 
We claim that 

, xeB{0,R), Ke 



l^'WI < (f 



A^-l 



s e N. (4.33) 



Indeed, if is: e and ii e e(-^) since i^Kue. e ii, we have 

\Ci\ = \{ni,-&KU£i) \ < W^KUliW <R- 



Next, using (4.32) and ||n,|| = 1, we have 



idriK) 



< 



|(n„jc)| + |c,-| ^ 2R 



(4.34) 



which readily implies ( |4.33 1. 

(B) We now use theorem [43 and corollary 48] to estimate the difference be- 
tween a Taylor polynomials and a Chung- Yao interpolation polynomial of a same 
function. To simplify, we omit the index s in the formulas. We have 

L[@;f]{x)-Ti-^{f){x) = [/(x) -Tj-^(/)(x)] - [f{x)-L[&,f]{x)] 

= E P^^\x)([0,...,0,x\nK,...,nK]f-[®K,x\nK,---,-tiK]f) 



(4.35) 



We call Si{x) and 52 (x) the terms in the above sum and prove that, for every 
X E B{0,R), both of them tends to as 5 — )■ oo. This will achieve the proof (since 
simple convergence on a compact set of nonempty interior implies convergence 
on any normed vector space topology on ^"^^^ (R^)). 
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(C) Since, in view of (4.33 1, the polynomials are bounded uniformly in s, 
that 5i (x) for X e b[o,R) follows from 



\0,..^.,0,x\ nK,..j ,nKjf- [&K,x\nK,--- ,^K]f 

d-N+l d-N+l 



^0 



(4.36) 



which is a consequence of the fact that the points of = 0*-^^ tend to together 
with the continuity of the divided differences of / as a function of the two groups 



of its arguments, see subsection 4.2 



(D) As for the term S2{x), since the right hand side goes to as * — )• oo, the 
conclusion follows from the following estimate. 



\S2ix)\< 



M 



{d~N+\)\ 



R 



1 + - 



|0||, xeB{0,R), (4.37) 



where ||0|| = ||0*^'^^|| := max{||tJ|| : ■& € 0| . To prove this, we observe that if 
N <i<d andKe (JJ'ij), the bound k.33) (in which M is replaced by e,_i) 
gives 



2R 

~5 



i-N 



xeB{0,R). 



(4.38) 



Moreover, for every a E B{0,R), using \\nK\\ < 1, we have 

/(^-^+i)(a)(x'^-',t?^u£,-,n7^)| <M||x||'^-'||t?^ui',-||||n^||'-^<Mi?'^-'.||0||. 
Hence, since vol(Aj_yv+i) = l/(t^ — 1)!, for x G B{0,R) we have 

>/[0,...,0,x] 

Combining the above estimates, we obtain 



< 



M 

{d-N+\)\ 



R'^^'\m 



(4.39) 



(4.40) 



^ f i — 1\/ 2i? 

I^2WI<£ 



i=N 



N-1 
M 



i-N 



M 



< 



{d 


-A^+1)! 




M 


{d 


-A^+1)! 




M 


{d 


-^+1)! 



{d-N+\)\ 



i-l 



N 



\®\\R 



,d-N 



i=N 

^d-1 



R'^^'\\@\ 



i-N 



5, 



10117?' 



,d-N 



1 + 



d-l 



(4.41) 
(4.42) 
(4.43) 
(4.44) 



This concludes the proof of theorem 3.1 
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4.5. An estimate on the error. The proof actually yields some estimate on the 
error between Chung-Yao interpolation polynomials and the Taylor polynomial 
at the origin. It is shown in the following corollary. 



Corollary 4.9. We assume that the assumptions of theorem \3.1\ are satisfied. If 
fe&-^+^{Q.) then 



max ||L[0i^■^/]W-T^--(/)W|| 

xtEB{0,R) 

= o(\\d^'^\\-i max ||/^"^+^)(a)||+ max \\f^'''^+^\a)\\\] . (4.45) 

V [aeB{0,R) aeB{0,R) } J 

where the constant involved in the symbol O does not depend on f. 

Proof. We turn to the term 5i (x) in the previous proof. For simplicity, we set 
m = d-N+ 1. Since / G C'"+^ (Q.), for all K E (^J , the mean value inequality 
gives 

I [0, ... , 0,x\nK, . . . , nK]f - [@K,x\nK, . . . , n/f]/| 

I I /^'"^ +1(0- ^) /^"^ (-^ + L ( Ok J - ) (n^) I 



< 



max 

A,„B(OM) 



Mrd^ < — max ||0||, (4.46) 

ml B{0,R) 



where &k = {Okj : z = l,...,m}. Using (4.37) and (4.46) in (4.35), we finally 
get 



max ||L[0,/](z)-T^-^(/)(x) 

xeB{0,R) 



< 



d 

N-l 



2R\ 



d-N+l 



i max \\f^'' 

(J-A^+l)!fi(o,i?)"^ 



+ 



1 



max 



{d - N + \)\ B{0,R 



-A'+l)|| j^d-N 



-N+2] 



d-l 



|0| 



1 + 



101 



= (Mi max \\f^''-^+^^\\+M2 max ||/^"^+^)||) ||0||. (4.47) 

V B{0,R) B{0,R) / 

□ 
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